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We consider an imperfect relativistic fluid, in which a shock wave develops and discuss its 
structure and thickness, taking into account the effects of viscosity and heat conduction in the 
form of sound absorption. The junction conditions and the nonlinear equations describing the 
evolution of the shock are derived with the corresponding Newtonian limit discussed in detail. 
As in the nonrelativistic regime, the thickness is inversely proportional to the discontinuity in the 
pressure. However, new terms of purely relativistic origin are also present. In particular, for a viscous 
polytropic gas, it is found that a purely viscous relativistic shock is thicker than its nonrelativistic 
counterpart, while for pure heat conduction, the contrary is true. 



I. INTRODUCTION 

The theory of relativistic shock waves was pioneered more than fifty years ago by Taub [Q, with the related 
junction conditions and adiabats further discussed by Israel |2j, Lichncrowicz |^ and Thorne These results were 
established for a relativistic perfect simple fluid, and since they do not involve any characteristic length scale, the 
shock front was described by a mathematical surface of zero thickness (abrupt transition). Many studies have been 
made extending these works to the nonlinear regime of relativistic hydrodynamics, as well as to ideal relativistic 
magnetohydrodynamics (e.g., see the book by A. M. Anile |^ and references therein). 

In this work, we are interested in the shock wave theory for an imperfect relativistic fluid. It is well known that 
for dissipative relativistic fluids, for scales smaller than the dissipation scale C associated to x/c, where x is the 
viscosity or heat conduction coefficient and c is the speed of light, ordinary Navier-Stokes formulae do not apply [1). 
We are concerned with fluid regimes whose characteristic lengths are larger than the dissipation scale, and therefore 
the classical theory for dissipative fluids may be used. 

All fluids are dissipative, and a nonrelativistic shock wave propagating in a dissipative medium cannot, in general, 
be considered an abrupt transition, but instead, as a region with a finite thickness. Its thickness is determined by 
the dissipative coefficients (i.e., viscosity and heat conductivity) In addition, shock waves propagating through 

a gas mixture that undergoes diffusion of one component, show similar characteristics as do those due to dissipation 

. This fact affects not only the evolution of the wave but can also have important effects on processes that depend 
upon the features of the shock wave. 

A pioneering study of this system was made many years ago by Koch ||l^ , which showed that if the shock velocity 
is greater than a given critical value, relativistic interaction of heat transfer and momentum transfer give rise to an 
increase in the velocity at the upstream end of the shock layer. The purpose of our work is to discuss several aspects 
of this system that were not considered in the work by Koch, thereby providing a more complete picture of relativistic 
shock waves. 

The main aim here is to derive the equation for the (generalized) Taub curve, as well as the general expression for 
the thickness of a plane shock wave in the weak relativistic regime, taking into account both the classical dissipative 
mechanisms (heat conduction, bulk and shear viscosity) and the associated sound absorption process. We shall not 
address the important issue of diffusion in a relativistic fluid or applications of our results to astrophysical processes, 
both of which will be left for future work. 

The paper is structured as follows. In the next section, we review briefly the Eckart formulation for a viscous and 
heat conducting relativistic simple fluid. In section III, we derive the junction conditions, as well as the corresponding 
expressions for the generalized Taub-Rankine-Hugoniot curves for a unidimensional shock wave in a nonequilibrium 
regime. The general expression for the thickness of the shock and its comparison with the nonrelativistic limit is 
derived in section IV . The Newtonian limits for the remaining expressions are discussed in the Appendix. The metric 
signature is (+,—,—,—). 
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II. IMPERFECT RELATIVISTIC FLUIDS: THE ECKART APPROACH 



The relativistic theory of imperfect fluids rests on two basic ideas. The first one is the local equilibrium hypothesis 
(LEH). It implies that for nonequilibrium fluids, state functions (such as entropy) depend locally on the same set of 
thermodynamic variables as do equilibrium fluids. In particular, the usual thermodynamic temperature and pressure 
concepts are maintained in the relativistic nonequilibrium regime. The second idea is the existence of a local entropy 
source strength (entropy variation per unit volume and unit time), which is always nonnegative, as required by the 
second law of thermodynamics. Mathematically, the LEH is represented by the Gibbs law, whereas the entropy 
law takes the form of a balance equation. Using these hypothesis in the fluid equations of motion, one finds an 
expression for the entropy source strength, as well as for the constitutive (phenomenological) relations. The perfect 
fluid equilibrium equations are recovered in the limit of a vanishing entropy production rate. However, an important 
point of difference in the treatment of relativistic and nonrelativisitic fluids by diff'erent authors should be stressed. 
In contrast to the Newtonian regime, in the relativistic domain, there exists an ambiguity related to the possible 
choices of the macroscopic hydrodynamic four- velocity. In the so-called Eckart formulation 12|, the four- velocity 



is directly related to the particle flux, while in Landau-Lifshitz's approach ||^, it is defined by the energy fiux. In 
principle, a general treatment should be able to deal with any of these "gauge" choices For simplicity and for 
the sake of a simpler comparison with previous studies, in what follows, we shall adopt the Eckart formulation. 

The thermodynamic state of a relativistic simple fluid is characterized by an energy-momentum tensor T"'^ , a 
particle current iV", and an entropy current S°' . The fundamental equations are expressed by the conservation laws 
(particles and energy-momentum) and the entropy flux equation, 

^.'^ = 0, rr = o, s>:^>o, (1) 

where iV^ is the particle flux, T'^'^ the stress tensor, and S"^ is the entropy flux (comma denotes space-time derivatives). 
In the Eckart frame, the particle flux and stress tensor can be written as 

iV^ = nu^', (2) 

T^"' ^ ^iu^'u" - ph^"" + TTh^"" + c-^q^'u'' + q^u^") + W , (3) 
with the entropy flux given by 



S^' = nkBCTu^' -Y- (4) 

The quantities n, p, p, a, T, and ks are the particle concentration, energy density, pressure, specific entropy (per 
particle), temperature, and the Boltzmann constant, respectively. The hydrodynamic 4-velocity is normalized 
according to it'^u^ = 1. The tensor 

/i^"^ = 5^" - u'^u'^ (5) 

is the usual projector onto the local rest space of u". The irreversible fluxes, tt, g^, and H^", are defined by 

^ = C0 , (6) 
q'' = Kh^^''{T,, -Ta,) , (7) 

n^"" = 77 +/i''"u^,a , (8) 

where k, C, and rj are the classical phenomenological coefficients (thermal conductivity, bulk and shear viscosity), and 
Oi/ = Uu^au"' is the four acceleration. The bulk viscosity stress, tt, represents an irreversible negative pressure, and 6 
is the scalar of expansion (divergence of 4-velocity). The heat flux q^^ is orthogonal to the 4-velocity, i.e., g^w^t = 0, 
whereas the shear- viscosity tensor H^", is symmetric, trace free, and space-like. 

For completeness, we recall that all dissipative fluxes tt, g'^, and H'^'^, as well as their space-time derivatives, are of 
first order of smallness in the equilibrium deviations. This is also true of the space-time derivatives of the reversible 
thermodynamic quantities n, p , p, a, and T. However, the source of entropy, i.e., the divergence of the entropy flux 

is a quantity of second order of smallness. In what follows, we write the 4-velocity as u'^ = 7(l,'i7/c), where 7 — 
(1 — w^/c^) is the Lorentz factor. 
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III. TAUB CURVES AND ENTROPY DENSITY CHANGE 



We now consider the junction conditions for a plane shock wave in a relativistic imperfect fluid and use them to 
derive the generalized Taub curve, as well as the associated entropy density change for weak shocks. The thickness of 
the shock wave due to the presence of viscosity and thermal conduction and taking into account the acoustic damping 
is, then, derived. 



A. Junction Conditions 



In an ideal fluid, the relativistic junction conditions are defined by the continuity equation for the particle current 
and the momentum and energy fluxes, i.e., the xx and Ox components of the energy-momentum tensor: 



[N^] =0, c [r°^] = 0, [T^^] = . 



(10) 



Square brackets denote the difference between the values of any of the mentioned quantities at large distances in 
front of the shock and inside it. We denote side 1 as the (far) upstream side. Choosing the spatial component of the 
four velocity along the a;-axis, it follows that u° = 7, = jv'^ /c. The non-null components of the projector tensor 
are: ft."" ~ — (u^) , /i"^ — —71*^ and h^^ — —7^. For convenience, the density particle current will be expressed as 
j = n^v. In this way, the conservation of the Oa; and xx components of the energy- momentum tensor take the form: 

(11) 
(12) 



n-f- = - = — 
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where we have introduced the specific enthalpy (per particle). 



(13) 



(14) 



and used u'q = u^dtu^ /c^, vP^ - 
uniform, i.e., all gradients vanish. 



u^dxU^ It has been assumed that at large distances from the shock, the flux is 



B. Generalized Taub curve 



In order to obtain the expression for the change in the entropy across the shock, we follow a procedure similar 
to that adopted by Thorne Q|. First we multiply (13) by (w/n+ Wi/ni) and then combine the result with p = 
n\ (uf ) (? = n^ [u^) c^, obtaining 



w^f + (p-pr) ( ^ + ^ ) ^ - {c+\v) 7^ ( - + ^ ) ^ 
' n ni I \ 3 1 \n ni I cn^ 



j dtn 
nc^ 7 



da: 
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j dxn 
n 7 



df.n 



(15) 



Multiplying (|12|) by (W7 + w 171), we get 
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(w^7^ - wljf) = - ( C + ) (w7 + wi7i) 



J 7 



j dfU 
771 



f 1 



1 c 
(w7 + W171) i — dtT + -jdxT 



T 



Finally, subtracting ([T^) from (|T^ and using 7^ = 1 

2 2 / 

W — W]^ = (p — pi 



cn J 
2 



n 7 



we obtain 



w ^ Wl 



+ I C + o'? ) 7^7iwi^- 



1 
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1 
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j d^n ' 
h Otn 

n 7 



dtn 



(16) 



(17) 



Equation (|l7|), together with the definition of j/c [Eq. (11)], are the generalized Taub junction conditions for a plane 
shock wave in an imperfect relativistic simple fluid. 



C. Weak shock wave: entropy density change 



We now consider the weak shock case, i.e., that for which all discontinuities are small. This means that differ- 
ences, such 'AS V — Vi, p — pi, etc., between the values in front of the transition layer and inside it are small. Thus 
differentiation with respect to x or ct increases the order of smallness by one, i.e., dV/dx is a quantity of second order 
of smallness. From (|lj), we see that the term involving the viscosity coefficients are of third order, while for heat 
conduction, the terms proportional to (w7 + W171) are of second order. The number particle conservation law can be 
written as 



n 



and considering that to lowest order, the temperature gradient satisfies p3| 



T 



dp^ 



the enthalpy density change, given by (^^, can be expressed as 
— = (p — pi) (w/n + vji/n) 



(18) 



(19) 



(20) 



+K (w7 + wi7i) <^ -dxT 
.3 



T 

7nc 



1- ^ 



dp 
9/i 



Note that the term proportional to T is proportional to c^^ through the dependence of 6 on the four velocity and 
time derivative (cf. eq. (p^). 

In the dissipationless regime, the resulting expression for the Taub adiabat is formally very similar to the Newtonian 
expression, as can be seen in Refs. |7| and For an imperfect relativistic fluid, the change in entropy in the transition 
layer is also of second order in the pressure, just as it is in the nonrelativistic case. However, as one may see from 
( po| ) three new purely relativistic terms come into play. In the Appendix, we show that eq. ( ^ ) yields the Newtonian 
expression previously found in the literature (e.g., Ref. |Q). 

To find the expression for the difference in the entropy density values far upstream and in the transition layer, we 
follow a standard procedure Q and develop w/n around its upstream value in powers of (p — Pi). We write the first 
law of thermodynamics as dvj~ dp/n + Tds, where s is the entropy per particle and then multiply by w, using the 
development of w/n. Keeping the zeroth order in wT in the second term and integrating, we get 



Wl 



2wiri(s-si) + 2— (p-pi) 



d_ 
dp 



ni 



(21) 



dp^ 
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As the derivatives of T and n are already of second order, we consider (w7 + W171) ~ 2 W171 in eq. (pO[). With this 
approximation and comparing with we obtain the entropy density change: 



Sl 



T 



(22) 



Therefore, as in the nonrelativistic case, the entropy density change is proportional to the heat conduction coefhcient. 
The nonrelativistic limit of this expression is trivial and coincides with the known expression |7| . 



IV. SHOCK WAVE THICKNESS 



Relativistic or nonrelativistic shocks are described by an evolving non linear wave. On the other hand, waves 
propagating in a viscous, heat conducting medium are damped. This fact can be phenomenologically described by 
an extra imaginary term in the dispersion relationship for the wave, i.e., by writing u ~ v^k — icLk^, where uj is the 
frequency, k the wavenumber, Vs the sound speed, and L is the absorption length (see Refs. 0,|l^). The equation we 
are seeking must be of the form 0| : 



d_ 

dt 



/ ~ Vsapf—f = cL—rf , 
ox I ox Ox'^ 



(23) 



where / is a suitable function that describes the wave profile. To obtain this equation, we shall follow a two-step 
procedure: we first find the nonlinear term (in the next subsection) and then proceed to find the quasi- acoustic 
damping contribution (in the subsequent subsection). 

A. Nonlinear term in shock waves 

In order to find the nonlinear contribution, we need only to consider equations for an ideal fluid. We consider now 
the local reference frame, in which the medium is at rest (comoving frame), and let 5v be a unidimensional velocity 
perturbation. We have (cf. Ref. [O) 



-irdv — dv 

& ot 



d_ 

dx 



Sv 



= 0, 



(24) 



dSv ^ d - 
^r— + dv—dv 
ot Ox 



Sv dp 
'^dt 



dx 



= 0, 



(25) 



Expanding /i and p, we have 



(26) 



p^Po + 6p, 



(27) 



where /zq and po are the background values. For a wave propagating to the left, we can write Sv = — [c? /v^^ Sp/ (/i -f p) 
and using d/dt = Vgd/dx in the second order terms, (pq) can be written as 



1 , , dSv d - 



-Sp—Sp 



(pto+Po) dx 



(28) 



and Eq. (24) as 



^pp+{^.o+Po)§^Sv = Vs 



{fio+Pa)v* \dp^ ) J dx 



(29) 



Deriving (pq) with respect to x and (E3) with respect to t and subtracting the resulting expressions, we get 
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where 



ai 



ifJ-o+Po) 



9V 



(31) 



a2 



(Mo +Po) 

Substituting d/dt = Vsd/dx and eliminating (9/9a; in both terms, we get 

Vs ot OX I ox 



where we define 



77 ^"1 



0.2 



in + p) vl 



(32) 



(33) 



(34) 



This equation has the same form as does the nonrelativistc one (see Ref. [Q ) , with (/i + p) replacing the rest mass 
density, the relativistic energy density derived twice with respect to the pressure, and the ratio of the sound speed to 
the light speed appearing explicitly. The last term in the square brackets is a purely relativistic correction. 

B. Dissipative term of the shock wave equation 

The acoustic relativistic damping length required by the complete nonlinear equation of a shock wave was derived in 
another context by Weinberg . We refer the interested reader to this work in order to see details of the derivation. 
Here, we just quote the final expression, 



cL = 



2{fi + p) 



1 1 



kT 



c2 a \dT 



(35) 



observing that Weinberg's expression is recovered when (Cp) — nc^ (ncp) and c = 1. The above expression is the 
same as the nonrelativistic result (see Ref. [Q), with {p + p) replacing the density of the rest mass and a relativistic 
correction proportional to kT appears explicitly. 



C. Solving for the thickness 



The complete equation of the evolution of a shock wave is obtained by adding to eq. (^3|) a term proportional to 
the second derivative with respect to x, which takes into account the dissipation. The final equation is then 



d a\. a dHp 

-TT-Vs-^ ]Qp- Vsapdp—dp = cL—r- 
Ot Ox I Ox Ox^ 



Following the usual analysis ||7|], we assume that Sp has the following dependence: 

6p^6p{^), i = x + Vwt 



where is the velocity of the wave. With this solution, Eq.(36) becomes 



d 



{vw - Vs) 6p - ^VsapSp"^ - cL^Sp 



0. 



(36) 



(37) 



(38) 



The solution to Eq.(||) is then 



P^ :^{Pi+P2) + :^{P2-Pi)tanh , 39 

where pi is the pressure far upstream and p2,the pressure far downstream. In the reference frame where the shock is 
at rest, we have for the pressure variation 

P - ^ (Pi +P2) = ^ (P2 -Pi)tanh , (40) 

where we have defined the "thickness" of the shock by 

A^— (41) 

Vsttp (P2 -Pl) 

We see that this expression is identical in form to the nonrelativistic one and proportional to the inverse of the pressure 
difference. The relativistic corrections are contained in the factors ap and cL. 



D. Analysis of the thickness 



In this subsection, we shall estimate the effect of the relativistic corrections to see if they increase or decrease the 
shock thickness. We shall examine them in the weak relativistic limit. We then write (see Appendix) 



and 



NR 



vIoiNR 



1 - 



A 



A 



NR 



V^aNR 



(42) 



(43) 



with Anh, X,aNR and rj given in the Appendix. The nonrelativistic expression for the shock thickness is 5 
Aa/ajsiR {p2 — Pl) where a — Ajvij/ff • We must evaluate 



A 

T 



cLa 



NR 



VsUpa vIoLNR Anr ' 

where the semi-equality holds for the weak relativistic case. Using the expressions in the Appendix, we find 

A ^ 2.?{l + W + 4K&)j 



= 1 



(44) 



(45) 
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"(C + |7?) {e+p) 




Ij 


pc^ 











For a more direct comparison of the relativistic thickness with the standard Newtonian result, we consider a polytropic 
gas and evaluate the above expression in two special cases: with viscosity alone and with thermal conduction alone. 



1. Polytropic gas 

In a classical polytropic gas, the energy density and enthalpy density are given by e = c^T = p/ {T — 1) and 
w = CpT = Tp/ (r — 1), where V = Cp/c^ = const, respectively. Hence {d'^e/dp^) ^ — and (dp/dT)^ = Cy (T ~ 1) 

and l/ci) ~ 1/cp = (r — 1)^ T/Tp. Replacing these formulae in the classical expressions for the internal energy s in 
eq. (^, we get 



Vs 27 


1 






I) 


c2 (r2 - 1) 


pc^ 









(46) 
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a. Only viscosity. If thermal conduction is absent, we find 

^ = 1 + ^^ (47) 

s c2(r + i) ^ ' 

We see that, in this case, the relativistic shock is thicker than the nonrelativistic one, with the increment proportional 
to the sound speed. 

b. Only heat conduction. If viscosity is absent, we obtain 

In this case, the relativistic shock is thinner than its nonrelativistic counterpart and the correction is again proportional 
to the sound speed. 



E. Entropy density change 



With the expression for the pressure given by (^^, we can express the entropy change as a function of the pressure 
discontinuity. We begin by writing explicitly the derivatives in the expression for 9, namely eq. (Uq) and replacing 
the time derivative with d/dt = Vgd/dx. Using dT/dx = (dT/dp)^ dp/dx + {dT/ds)pds/dx ~ {dT/ dp) ^ dp / dx and 
dn/dx — {dn/dp)^dp/dx + {dn/ds)pds/dx ~ (dn / dp) ^ dp / dx in Eq.(|2^) and evaluating dp/dx from (|40|), we obtain 
the following expression for the entropy density change in a reference system in which the shock is at rest: 



K 



71 

P {P2 



] 



dT 
dp 



^n^c^ 



1 



dp 



dn 
dp 



ScL cosh^ (x/A) 



(49) 



where the factor c in cL docs not add an extra power in the speed of light (see Appendix) . In the nonrelativistic case 
, the entropy reaches a maximum inside the shock |^ and is of second order in the pressure discontinuity. 



V. CONCLUSION 



In this paper, we have extended previous studies of shock waves done in the nonrelativistic domain to the weak 
relativistic case. Considering dissipative relativistic fluids in the range of validity of the Navier-Stokes-Fourier theory 
we have obtained expressions for the entropy density change and the shock thickness that coincide in form with 
nonrelativistic ones. In each of the factors in the equations, purely relativistic corrections appear explicitly. We studied 
the expression for the shock thickness for a polytropic gas and analyzed the effect of corrections in two important 
limits defined by the presence of viscosity or heat conduction. When only heat conduction is taken into account, 
the relativistic shock is thinner than for the nonrelativistic case. This result can be understood by observing that 
heat conducting fluids can develop "thermal discontinuities" [QJs), i.e., they allow for discontinuities in the velocity, 
pressure and density of the fluid flow, while the temperature remains constant. On the other hand, when only 
viscosity is present, the shock thickness is larger than for its nonrelativistic counterpart and, hence, the tendency to 
erase singularities is stronger in the relativistic limit than in the Newtonian one. This difference in the effect of the 
relativistic corrections can also be understand, as follows. Viscosity provides the mechanism to convert a portion of 
kinetic energy of the gas flowing into the discontinuity into heat. This conversion is equivalent to the transformation of 
the energy of ordered motion of gas molecules into energy of random motion by the dissipation of molecular motion. In 
this respect, heat conduction has an indirect effect on the conversion process since it only participates in the transfer 
of the energy of random motion of the molecules from one point to another, but does not directly affect the ordered 
motion. The corresponding relativistic corrections seem to amplify these effects. 

A comment on the entropy change is in order. When the pre-shock gas has a low temperature, we are in the 
strictly Newtonian limit and, in this sense, the fact that the entropy density change reduces to its Newtonian analog 
is equivalent to requiring that the theory has the correct low-speed limit. This fact contains no new information. But 
when the pre-shock fluid has relativistic internal speeds, the shock weakness does not imply a Newtonian propagation 
velocity of the shock and, hence, this case is not covered by the Newtonian treatment. In this sense, the result that we 
have obtained, that the entropy change still reduces to the Newtonian expression is new and potentially interesting. 
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Finally, it should be mentioned that although first-order theories are successful in revealing the physics underlying a 
large class of phenomena, they present some experimental and theoretical drawbacks. In its classical version, the linear 
constitutive equations (||)-(|8|) are not adequate at high frequencies or short wave lengths, as manifested in experiments 
on ultrasound propagation in rarefied gases and on neutron scattering in liquids |15| . In additon, they also allow for 
the propagation of perturbations with arbitrarily high speeds, which although unsatisfactory on classical grounds, 
is completely unacceptable from a relativistic point of view. Furthermore, they do not have a well-posed Cauchy 
problem and their equilibrium states are not stable. Several authors have formulated relativistic second-order theories 
which circumvent these defficiencies |l6| , pT| , p^|Jl9| ] . In a forthcoming paper, we intend to extend our considerations to 
this class of theories. 
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VI. APPENDIX 



In this appendix, we obtain the nonrelativistic limits of the magnitudes discussed in this paper. We begin with 
the nonrelativistic limit of 7^(/i-|-p), neglecting the pressure since nmc^ 3> p in the nonrelativistic limit. Thus, 
7^ (/i + p) fJ. = ^"^nrac^ -f 7^£, where £ is the internal energy density, i.e., the energy associated with internal 
degrees of freedom. In the limit of small velocities, nm — > ^, where p is the mass density and, therefore, ^'^nmc? 

pf? -f- (l/2)pu^. Taking 7 = 1 in the expression for the internal energy density, we 



7pc^ = (l — u^/c^) ^^"^ p(? 



obtain the desired limit: 7^ {p. + p) pc^ + (l/2)pw^ -I- pe, where e is the internal energy per particle. 



A. Nonrelativistic limit of the Taub curve 



At the limit c — > 00 in Eq.(EG), the terms in the square brackets can be neglected. Thus, 



w Wi 



w - — {p - pi) ! + - (w + wi) d^T 



n n\ ) J 



(50) 



w — [p — pi) tWOxT = Wi + [p — pi) h— WiO^T. 



n J 



m J 



(51) 



Taking the square-root, 



l-ip-pi) — -—dxT 
wn jw 



1/2 



= Wl 



1 + {p ~ Pi) ^— + —d^T 
Wini jwi 



1/2 



(52) 



Assuming that {p — pi) and dxT are small, we have 



w 1 - 7: (p-Pi) 



1 1 K 
wn 2 jw 



dxT] =wi 1 + -(p-pi) 



1 1 K 



wini 2jwi 



dxT 



(53) 



or, rearranging terms. 



1 , , /I 1 \ 

w - Wl = - (p -pi) - H + -d^T, 

2 \n nij J 



(54) 



which is the standard nonrelativistic expression for the Taub Curve 
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B. Nonrelativistic limit of the shock thiclcness 



The nonrelativistic limit of the shock thickness is derived from Eqs (34) and (35), the expressions for nonrelativistic 
ttp and L, respectively. Using the expression for 7^ + p), derived in the introduction to the Appendix, we rewrite 
eq. (p3) for ap as 



1 



5^ / ^ £ \ 



vl {pc^ + e + p) 



(55) 



In the weak relativistic limit, we obtain 



oip = aNR - r] , 



(56) 



with 



,,2 r 



aNR 



(d^ 

pvf V 9p'^ 



vlotMR , 



(57) 



Pc I pvj 



2 ^ \dp^ 



(58) 



It is convenient to express the dispersion relationship as fc = jw/vs + ij'^cLuj'^/v'^, where L is defined in Eq.(35). 
Using {p+p) /(? = p + [e + p) /c^, in the weak relativistic hmit, we have 



cL — Apfj^ — A = A 



NR 



A 



A 



NR 



(59) 



where 



^NR - TT 

2p 



K / 1 1 



2 V c. 



(60) 



A = 



2p2c2 



^ 4 \ . , kT 2 { dp 



(61) 



with Cy — Cy/p. the specific heat per unit mass. 

Using Eqs.(|5q)-(|6l|) in eq. (^) and taking the limit c ^ 00, we obtain the standard expression for the thickness 
of a nonrelativistic shock lf7|: 



5 = 



(92y/V)Jp2_pi) 



(62) 
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